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We investigate meson mixing and meson coupling constants in pion superfluid in the framework 
of two flavor NJL model at finite isospin density. The mixing strength develops fast with increasing 
isospin chemical potential, and the coupling constants in normal phase and in the pion superfluid 
phase behave very differently. 
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I. INTRODUCTION 

Recently, the study on the phase diagram of Quantum 
Chromodynamics (QCD) is extended from finite tem- 
perature and baryon chemical potential to finite isospin 
chemical potential. The physical motivation to study 
pion superfluid formed at high isospin density is related 
to the investigation of neutron stars, isospin asymmetric 
nuclear matter and heavy ion collisions at intermediate 
energy. 

Different approaches, such as chiral perturbation 
theory [lL El, M : Eh Nambu-Jona-Lasinio (NJL) modelQI 
HI, B, M IsL lid. random matrix methodfl^]. lattice 
QCDQJ, ladder QCD [3] and strong coupling lattice 
QCD '15j, have been used to investigate the QCD phase 
structure at finite isospin density. It is found that there 
is a phase transition from the normal phase to the pion 
superfluid at a critical isospin chemical potential which 
is exactly the pion mass in the vacuum, fi c j = m n . The 
temperature and chemical potential dependence of me- 
son masses in the pion superfluid is calculated in the NJL 
model[|H, linear sigma model [l6|, [TtJ and <fi 4 theorv[l8j. 
The familiar meson properties are significantly changed 
in the pion superfluid phase. In the normal phase which 
is controlled by spontaneous chiral symmetry breaking 
and its restoration, the scalar and pseudoscalar collective 
modes of a two flavor quark system are a and tt mesons. 
However, they are no longer the eigenmodes of the Hamil- 
tonian in the pion superfluid, and the new eigenmodes 
are the linear combination of them. It is the combina- 
tion that guarantees the Goldstone mode corresponding 
to the spontaneous isospin symmetry breaking Q. 

To investigate the dynamic processes which may be 
considered as signatures of the pion superfluid, one needs 
not only the meson masses but also other in-medium me- 
son properties in the superfluid. In this Letter, we will 
study the hadron mixing and coupling constants for a 
and 7T mesons at finite isospin density. 

It is well known that, the perturbation theory of QCD 
can well describe the new phases at extremely high tem- 
perature and high density, but the study on the phase 
structure at moderate temperature and density depends 
on lattice QCD calculation and effective models with 
QCD symmetries. One of the models that enables us 
to see directly how the dynamic mechanisms of chiral 



symm etry breaking and restoration operate is the NJL 
model[l9[ applied to quarks [2p|. When the parameters 
in the model are properly chosen, the calculated chiral 
phase transition line [2d. l21j in temperature and baryon 
chemical potential plane is very close to that obtained 
with lattice QCD[22|]. Recently, this model is used to in- 
vestigate color superconductivityj23j at moderate baryon 
density and pion superfluidity [EIlJE II, B El EH at finite 
isospin density. We will study the meson mixing and cou- 
pling constants in the pion superfluid in the framework 
of this model. 

The Letter is organized as follows. In Section [TT] 
we review the NJL model in mean field approximation 
for quarks and random phase approximation (RPA) for 
mesons in the normal phase at finite isospin chemical po- 
tential. In Section InTl we focus on the meson mixing and 
meson couplings to quarks in the pion superfluid. We 
summarize in Section llVI 



II. MESONS IN NORMAL PHASE 

The mesons in the normal phase at finite isospin chem- 
ical potential are controlled by chiral dynamics and ex- 
plicit isospin symmetry breaking. In this section we re- 
view the quark propagator in mean field approximation 
and meson polarizations in RPA in the NJL model at 
finite temperature and baryon and isospin chemical po- 
tentials. The Lagrangian density of the two flavor NJL 
model at quark level is defined as(20j 

C = 4) (z7 M <9 M - to + fijo) ijj + G (V'V') 2 + (4>ij5nip) 

with scalar and pseudoscalar interactions correspond- 
ing to a and n excitations, where too is the cur- 
rent quark mass, G is the four-quark coupling con- 
stant with dimension (GeV)~ 2 , (i = 1,2,3) are the 
Pauli matrices in flavor space, and /i = diag (p u , fid) = 
diag (fiFj/3 + fJ-i /2, iib/3 — fJ-i /2) is the quark chemical 
potential matrix with fi u and fid being the u- and <i-quark 
chemical potentials and fiB and fii being the baryon and 
isospin chemical potentials. 

At zero isospin chemical potential, the Lagrangian den- 
sity has the symmetry of U B (l) ® SU I (2) (g) SU A (2) cor- 
responding to baryon number symmetry, isospin sym- 
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metry and chiral symmetry. At finite isospin chemical 
potential, the isospin symmetry SUi(2) and chiral sym- 
metry 5C/a(2) are explicitly broken to Ui(l) and Ua{1), 
respectively. Therefore, the chiral symmetry restoration 
at finite isospin chemical potential means only degener- 
acy of a and ttq mesons, the charged 7r+ and 7r_ behave 
still differently. 

Introducing the chiral and pion condensates 

a = (ipip) = a u + (T d , 

a u = (uu), a d = (dd), 

7r = ^(-tpi^r+ip) = V2(^i75r_^) (2) 

with t± — (ri ± «r 2 )/v / 2 and taking all the condensates 
to be real, the quark propagator in mean field approxi- 
mation can be expressed as a matrix in the flavor space 



S(p) = 



with the four elements! 
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where = E ± =p /is/3 arc the energies of four 



quasi-particles with E ± = <J (E ± //j/2) + 4G 2 tt 2 , £ 



p| 2 + Mq and effective quark mass M q = mo — 2Gr, 

and A± = [1 ± 70 (7 • p + M q ) /E] /2 are the energy pro- 
jectors. The quark propagator ([3]) is the background of 
our following calculations for quarks and mesons. From 
the definitions of the chiral and pion condensates © , it 
is easy to express them in terms of the matrix elements 
of the quark propagator, 



p 



\Tr [iS uu (p)] , 

o"d = - ^2 Tr [iS d d(p)} , 
p 

7r = ^2Tr[(S ud (p) +S du (p))^ 5 ] , 
p 



(5) 



where the trace Tr = TrcTrr) is taken in color and 
Dirac spaces, and the four momentum integration is de- 
fined as J2 P = iT J2 n J d 3 p/(27r) 3 with p = iu n = 
i(2n + 1)ttT (n = 0, ±1,±2, •••) at finite temperature 
T. Note that, for [Lb — or = there is always 
Cu = o'di since in this case the chemical potential differ- 
ence between u and u is the same as that between d and 



d. The coupled set of gap equations ([5]) determines self- 
consistently the three condensates. The isospin chemical 
potential dependence of the effective quark mass M q and 
pion condensate 7r is explicitly shown in The phase 
transition from the normal phase to the pion superfluid 
is of second order. 

In the NJL model, the meson modes are regarded as 
quantum fluctuations above the mean field. The two 
quark scattering via meson exchange can be effectively 
expressed at quark level in terms of quark bubble summa- 
tion in RPA[20(. The quark bubbles, namely the meson 
polarization functions are defined as [2 11] 

Tlmn(k) =iJ2 T r (r* m S(P + k)T n S{ P )) (6) 



with the trace Tr = TrcTrpTrT) taken in color, flavor 
and Dirac spaces and the meson vertexes 



m = a 



r — 

1 n) — 



I 

ir + 75 m = 7T+ 

ir_75 m = 7r_ 

?r 3 75 m = ttq . 



m = a 



ri = 



1 

IT_75 m = 7T + 

ir + 75 m = 7r_ 

ir 3 7 5 m = n Q 



The explicit T, fis and /i/ dependence of the meson po- 
larization functions ([B]) used in the following discussion 
can be found in Appendix B of Q. 

At fij < /ij where [i c j is the critical isospin chemical 
potential for the pion condensate, the system is in the 
normal phase with diagonal quark propagator, and the 
bubble summation in the effective interaction in RPA se- 
lects its specific isospin channel by choosing at each stage 
the same proper polarization[20]. Therefore, the meson 
masses M m (m = a, 7r + , 7r_, no) which are determined by 
poles of the meson propagators at fco = M m and k = 
are related only to their own polarization functions [20( , 



1 — 2GII mm (fco)|fc =M ro — 0- 



(8) 



From the comparison of these mass equations with the 
gap equation for the pion condensate tt at T = [Ib = 
but finite fj,i, the critical isospin chemical potential where 
the normal phase ends and the pion superfluid starts is 
exactly the pion mass in the vacuum M, [l3| . /if = 77U-, 
and the /zj dependence of the meson masses is simple 
M„(m) = m a , M Wo (m) = and M w± (/ij) = T 
fix, where m a and are the a and 7r masses at T = 
Ms = Mi = 0. The isospin neutral mesons keep their 
vacuum masses, while the isospin charged mesons change 
their masses linearly. These relations hold until the pion 
condensate starts. 

The meson couplings to quarks, g m qq, are related to 
the residues at the corresponding poles of the meson 
propagators (20j. 
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[dn mm (k )/dk 2 ] 



=M„ 



(9) 
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III. MESONS IN PION SUPERFLUID 

In the phase with finite pion condensate which leads to 
spontaneous isospin symmetry breaking, the quark prop- 
agator contains off-diagonal elements, we must consider 
all possible isospin channels in the bubble summation 
in RPA[20|. While there is no mixing between ttq and 
other mesons U„ oa (k) = U 7T(>7T+ (k) = U 7T(>7T _(k) = 0, 
the other three mesons are coupled together, and the ef- 
fective interaction in two quark scattering via exchanging 
these mesons in RPA becomes a summation in the meson 
space, 

U(k) = T* m M mn (k)T n , to, n = cr, 7r + , 7r_ (10) 
with the meson matrix M(k) defined by 



M(k) = 



2G 



Of 

" T M(k), (11) 



1 - 2GTl(k) D(k) 
where 1 — 2GIT(fc) is the meson polarization matrix 
1 - 2GU 

1 — 2GTl aa 2GII (T7r+ — 2GIT ( j ir 
— 2GII 7r+(T 1 — 2GII Tr+7r+ — 2GH 7V+7r _ 
— 2GTL 7 r_(j — 2GiT 7r7 i-^ 1 2GT\-^ _^ _ 

D{k) is its determinant, 

D(k) = det (1 - 2Gn(fc)) , 



(12) 



(13) 



and M(k) is defined as M(k) = D(k)/(1 - 2Gn(fc)). 
In this case, cr, 7r + and 7r_ are no longer the eigenmodes 
of the Hamiltonian, the new eigenmodes are linear com- 
binations of them. In the following we call these new 
eigenmodes in the pion superfluid phase as a, 7r + and 

7T_. 

The ttq mass and coupling constant are still controlled 
by its own polarization function at k = 0, 

1 — 2Gn 7ro7ro (fco)|fe 0= Af 7ro = 0, 

<4 g _ = [dU^ (k )/dk 2 ] _1 | feo=M . , (14) 

since it is independent of the other mesons. At T = [1b = 
0, the 7To mass is exactly equal to the isospin chemical 
potentialQ, M Wo (fj,/) — /ij. At the critical point pf\ 
;u - . it is continuous with the solution M 7TQ — in the 
normal phase. 

The masses of the new eigenmodes a, 7F4. and 7r_ are 
determined by the pole of the effective interaction at k = 
0, 



D(k )\ ko= M„ =0, 9 = a, 7T+, 



(15) 



It can be proven analytically that there is always a zero 
solution which guarantees the Goldstone mode, M^, = 
0, corresponding to the spontaneous isospin symmetry 
breaking [8(. 



In order to derive the coupling constants for the new 
modes, we first expand the effective interaction U around 
the meson mass fcg = Mg at k = 0, 



U(k 



_2G r^ mn (M e )r n 

2M. 1.2 _ A/ r2 > \ W ) 



{dD{k a )/dkl)\ kQ=Me kl-Ml 



and then make a transformation from the coupled me- 
son space (cr, 7r + , 7r_ ) to the diagonalized meson space 
(<7, 7f + , 7f_ ) . With the help of the pole equation ([Tij]) for 
the mass Mg, we can derive the relations between the 
diagonal and off-diagonal elements of the matrix M. , 



MmmMnn - M mn = (1 - 2GII„).D = 0, 

I, rn, n = a, 7r + , 7r_, I m =/= n, (17) 



where all the quantities are evaluated at ko = Mg. Tak- 
ing into account the symmetric property of M, 



M ■mil 



(18) 



the effective interaction U via changing a 0-meson in the 
diagonalized meson space can be written as 



U(k ) 



2GM(M 9 ) 



r;r fl 



[dD{k Q )/dkl)\ ko =M B k 2 - Ml 



(19) 



in terms of the new meson vertex 

r 



M I n in (M e ) T m JM(M 9 ) 



r *e - E \/M mm (M ) T* m /^M(M ), (20) 

m 

where M(M e ) is defined as M(Mg) = T, m M mm (M e ). 

Since the meson coupling constant is defined as the 
residue of the effective interaction at the pole, we can 
extract the 0-meson coupling constant gg q g from (fH 
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9 6 qq 



2GM (Mg 



" (dD(k Q )/dkl)\ ko=Me 



(21) 



To explicitly describe the meson mixing in the pion 
superfluid phase, we can introduce mixing angles between 
two mesons, for instance, the angels a between 7r + and 
cr, (3 between 7r_ and a and 7 between 7r+ and 7r_ in the 
cr-meson channel, 



tana = \ M 7I+7I+ (M W ) / J M aiJ (M w ) 



tan/3 = J M„_„_(M W ) / d Maa(M w ) 



tan 7 = yj M„ + ^ + (M w ) j M„_ % _(M W ) 

= tana/ tan (3. (22) 

It is easy to see that only a and (3 are independent. The 
mixing angels in the 7f + and W _ channels can be defined 
in a similar way. 
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FIG. 1: The phase diagram in \ii — fis plane at fixed temper- 
ature T. 



Since the NJL model is non-renormalizable, we can 
employ a hard three momentum cutoff A to regular- 
ize the gap equations for quarks and pole equations for 
mesons. In the following numerical calculations, we take 
the current quark mass mo = 5 MeV, the coupling con- 
stant G = 4.93 GeV~ 2 and the cutoff A = 653 MeV [p. 
This group of parameters corresponds to the pion mass 
TOjt = 134 MeV, the pion decay constant f n = 93 MeV 
and the effective quark mass M q = 310 MeV in the vac- 
uum. 

The phase diagram in the fij — fiB plane at fixed tem- 
perature T is shown in FigflJ The system is in the normal 
phase at low /ij or high /i b and in the pion superfluid at 
high fii and low fis- At zero temperature, the critical 
isospin chemical potential fi c j does not change with [Ib 
till 



Hb = 3(M q - m^/2) - 730 MeV 



(23) 



which is determined by the comparison of the gap equa- 
tion for the pion condensate at ir = and T = with 
the pion mass equation in the vacuum. With increasing 
temperature, the critical value n\ increases and the pion 
superfluid region decreases, due to the melting of pion 
condensate in hot mediums. 

The mixing angles a, (3 and 7 in the cr-meson channel 
are shown in the upper panel of Fig[2] as functions of \xi 
at T = jitg = 0. In the normal phase with /i/ < /ij, a 
and 7r themselves are the collective excitation modes of 
the system, and there is no mixing among them. In the 
pion superfluid phase with /i/ > /Ltf, a and (3 indicate the 
7T_|_ — (7 and 7r_ — a mixing strengths in the cr-meson chan- 
nel, and 7 reflects the relative strength between them. 
While in the very beginning of the superfluid any mix- 
ing is weak, it develops fast. For /ij > 200 MeV which 
corresponds to a = 7r/4, the mixing is already so strong 
that the contribution from 7r+ to a is larger than that 
from a. Similarly, the 7r_ -component in a becomes more 
important than the cr-componcnt itself even for /i/ > 150 
MeV. Therefore, at not very high isospin chemical poten- 
tial the 7r_- and 7r + -components start to dominate the a 
mesons. While the 7r_ — a mixing is always stronger than 
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FIG. 2: The mixing angles a, /3 and 7 in the <r-meson channel 
as functions of /ir at T = /xg =0 (upper panel) and T = 
O.lGeV and (ib = 0.4GeV (lower panel). 



the 7r + — a mixing, namely /3 > a, the relative strength 
shown by 7 decreases with increasing isospin density. 

At zero temperature, since the phase structure of the 
pion superfluid does not change till the baryon chemi- 
cal potential \ib = 730 MeV, see the vertical straight 
line in FigfTJ the mixing angles at finite baryon chemi- 
cal potential /is < 730 MeV are exactly the same as the 
ones shown in the upper panel of Fig[2j Since the pion 
superfluid happens only at low temperature, we calcu- 
lated the fii dependence of the meson mixing angels at 
T = 0.1 GeV and /is = 0.4 GeV, shown in the lower 
panel of Fig|2] From the comparison with the upper 
panel, the shape of the mixing angles at finite temper- 
ature and baryon chemical potential is almost the same 
as that at T = \ib = 0, the only remarkable change 
is the critical value /ij. It increases from 134 MeV at 
T = ii B = to 153 MeV at T = 0.1 GeV and t i B = 0.4 
GeV. 

From our numerical calculation, the mixing angles in 
the ?F_ -meson channel behave similarly. However, the 
case is significantly changed for the Goldstone mode W+ . 
The vertex for W+ can be greatly simplified as 



iv + = (TV - r 7r _) /v2 



(24) 



at any temperature and baryon and isospin chemical po- 
tentials. The Goldstone mode contains only tt + - and rr_- 
components, and the fractions for the two components 
are exactly the same. 
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FIG. 3: The coupling constants for cr, no, 7T_ in the nor- 
mal phase and a, 7ro,7f+, 7f_ in the pion superfluid phase as 
functions of fii at T = /j,b = 0. 



exists no a meson, and the coupling constant g^qq drops 
down to zero at the critical point /ij and keeps zero at 
Hi > fij. For the other three mesons 7To,7f+ and 7f_, their 
masses are below the threshold value, and the coupling 
constants are nonzero. 

From Fig [21 the meson mixing angels are only slightly 
changed by finite temperature T and baryon chemical 
potential /is, when T and \ib are in a reasonable region. 
Therefore, it can be expected that the meson coupling 
constants will not be significantly changed by finite T 
and fiB , except a remarkable shift of the critical isospin 
chemical potential fif from 134 MeV to a higher value. 



IV. SUMMARY 



The meson couplings to quarks are shown in Fig|3] as 
functions of fj,i at T = /zb = 0. In the normal phase 
with [ij < /ij, the coupling constants are calculated 
through the diagonal polarization functions II mm . At 
T = ii b = 0, n mm depends only on the quark mass M q , 
meson mass M m and isospin chemical potential /ij. Since 
M m for isospin neutral mesons a and ttq and M q are con- 
stants, the coupling constants g a qq and g^ oq q are [Li inde- 
pendent. In the pion superfluid phase, the coupling con- 
stants, determined by the diagonal and off-diagonal po- 
larization functions n m „, behave very differently. From 
the mass relation — > M no at \i\ — » oo, the couplings 
9n-qq an d g-n Q qq approach each other at large enough /!/. 
Note that g^ Q qq is no longer a constant but changes slowly 
in the pion superfluid. Since we did not consider the me- 
son widths in the pole equation (fT5j) , the condition for 
a meson to decay into a q and a q is that its mass lies 
above the q — q threshold. The meson masses, calculated 
through the pole equations ([5]) in the normal phase and 
(fT4t and (fl~5| in the pion superfluid, are shown in Fig. 15 
of [8| . In the pion superfluid, the heaviest mode is a and 
its mass is beyond the threshold value. As a result, there 



We have investigated the meson mixing and meson 
coupling constants in the NJL model at finite isospin 
chemical potential. In the pion superfluid phase, the 
normal mesons are no longer the collective excitation 
modes of the system, and the mixing among them be- 
comes important. For the Goldstone mode, it contains 
only charged pions and their fractions are exactly the 
same in the whole superfluid region. For the other new 
eigenmodes, the meson mixing starts to control the sys- 
tem at /ij > 150 MeV which is only a little bit higher than 
the critical value ^ = m n = 134 MeV at T = fi = 0. 
The coupling constants for the conventional mesons in 
the normal phase and for the new eigenmodes in the pion 
superfluid phase behave very differently. The splitting of 
meson mass and coupling constant due to explicit and 
spontaneous isospin symmetry breaking can be used to 
further calculate the itit scattering and the number ratio 
7r + /7r_ in isospin asymmetric matter, which can help us 
to measure the meson properties at finite isospin density. 
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